Pattern forming instability induced by light in pure and dye-doped nematic liquid 

crystals 
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We study theoretically the instabilities induced by a linearly polarized ordinary light wave incident 
at a small oblique angle on a thin layer of homeotropically oriented nematic liquid crystal with 
special emphasis on the dye-doped case. The spatially periodic Hopf bifurcation that occurs as the 
secondary instability after the stationary Freedericksz transition is analyzed. 



I. INTRODUCTION 

Liquid crystals (LCs) are known to demonstrate a very 
rich variety of interesting optical phenomena which have 
been studied intensively during the last two decades. A 
nematic LC behaves optically as a uniaxial anisotropic 
medium with the optical axis along the local molecular 
orientation described by the director n(r,t). Moreover, 
when light propagates through the nematic, its electric 
field exerts a torque on the molecules which can induce 
molecular reorientation. When there is only an ordinary 
wave in the LC (polarization perpendicular to the plane 
containing the optical axis and wavevector), the initial 
distribution of the director becomes unstable when the 
intensity of light reaches a certain critical value. This 
is the so-called Light Induced Freedericksz Transition 
(LIFT). The director reorientation leads to a change of 
birefringence and, as a consequence, the polarization is 
changed as light propagates through the layer jjl || . 

It is known ^ ^] that a periodic equilibrium configu- 
ration of the nematic director can appear in a thin-film 
LC in the magnetic or electric field induced Freedericksz 
transition under certain conditions. Our work is devoted 
to the search of analogous phenomena in the LIFT of ne- 
matic LCs including the dye-doped case. We will show 
that the Hopf bifurcation that occurs as a secondary bi- 
furcation after the LIFT leads indeed to a periodic pat- 
tern, although the mechanism is here quite different (see 
Conclusions). Doping is important because the LIFT 
threshold of a dye-doped nematic can be two orders of 
magnitude smaller than for a pure nematic. The nature 
of this anomalously low threshold was the subject of nu- 
merous studies (see ||, ^| and references therein). The 
fact that the threshold intensity is low allows the spot 
size of the light to be much larger than the thickness of 
the layer, thus a large aspect ratio system can be realized. 

This paper is organized as follows. In Sec. II we 
present the theoretical description of our problem. In 
Sec. [II we perform the linear stability analysis of 



the homeotropic state which gives the threshold for the 



LIFT. The numerical method of calculating the station- 
ary distorted state is described in Sec. IV. Finally, in 
Sec. |v|, we do the linear stability analysis of the station- 
ary distorted state with respect to general perturbations 
in the plane of the nematic layer. 



II. THEORETICAL MODEL 

We consider a linearly-polarized plane wave incident 
at a small oblique angle /3o on a layer of dye-doped ne- 
matic LC which has initially homeotropic alignment (see 
Fig. 1). The light is polarized along the y-axis i.e., 
we deal with an ordinary wave. Strong anchoring of 
the nematic at the boundaries of the layer is assumed. 
The optical torque acting on the director is given by 
t = £ e ///167r(n • E*)(n x E) + c.c, where E is the 
amplitude of the optical electric field, £ e // = e a + (■ 
e a = £ii — e± is the dielectric anisotropy and e± (e\\) is 
the dielectric permittivity (at optical frequency) perpen- 
dicular (parallel) to n. ( phenomenologically describes 
the effect of certain dye dopants (£ e // = £a in a pure 
LC) and can be both positive and negative and depends 
on dye concentration, molecular structures of both host 
and dye materials, on the wavelength of light, and on 
the temperature || ^J. Obviously, the electrical part of 
the free energy will contain the same factor £ e //- The 
density of the free energy of the dye-doped nematic LC 
is thus assumed in the form: 



F = F, 
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where F elastlc = K x /2 (V • n) 2 + K 2 /2 (n • V x n) 2 + 
K 3 /2 (n x V x n) 2 is the standard Frank free energy den- 
sity and Ki , K 2 , K 3 are respectively the splay, twist and 
bend elastic constants of the LC • 

We first assume that the director components depend 
only on z,t. We introduce the angles 0(z,t) and (p(z,t) 
(see Fig.l) so that n= (sin 8, cos 9 sin <p, cos 9 cos <p). Us- 
ing the standard variational principle JD and taking the 
dissipation function in the form R = -^n 2 , where 7 is an 
effective rotational viscosity, the equations of motion for 
9{z,t) and ip(z,t) can be derived. 

In addition we need Maxwell's equations to determine 
the electric field. These equations contain the complex 
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FIG. 1: Geometry of the setup: linearly polarized light along 
the y-direction incident at angle /3o on a nematic LC layer 
with the director no || z (homeotropic state). The compo- 
nents of the director n are described in terms of the angles 
8,ip (9 — if — in the homeotropic state). 



dielectric tensor which depends on the director compo- 
nents: 



(2) 



In (g) 7a = 7± — 7||, where 7^ and 7|| are the imagi- 
nary parts of the dielectric permittivity for E perpen- 
dicular and parallel to n respectively. They describe 
the absorption effect by the dye, so they vanish in pure 
LCs. The magnetic anisotropy at optical frequencies 
can be neglected. Since the components of the dielec- 
tric tensor depend on the z-coordinate only, we may 
use the stratified medium approach for describing wave 
propagation We write the electric and magnetic 

fields in the form: E(r,t) = l/2(E(z, t)e is ° k ° x e- iut + 
c.c), H(r,<) = l/2(U(z,t)e isokoX e- luJt +c.c), where k = 
u)/c is the wavenumber in vacuum and sq = sin(/3o). Here 
E(z, t), H(z, t) are amplitudes that vary slowly in time 
compared to lu^ 1 and obey the equation: 



d<3> 

— = ifc D*, 

dz 



where 
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( 5 ) 

The z-component of the electric field can be found from 
the following relation: 



E z 



so 



-E x 



-yz 



(6) 



We will examine the case £ e // > so that the preferred 
orientation corresponds to the director parallel to the 



electric field n||E. Since in our geometry initially n _L E, 
the homeotropic state will cease to be stable above some 
critical intensity of the incident light. The reorientation 
of the LC leads to modification of the electric field polar- 
ization inside the LC owing to the fact that it becomes 
an inhomogeneous anisotropic medium. 



III. STABILITY ANALYSIS OF THE 
HOMEOTROPIC STATE 

We first perform the linear stability analysis of the 
homeotropic state (8 = ip = 0). The linearized equation 
of motion for ip(z, t) has the following simple form: 



jd t ip = K 3 d z ip ■ 
E\ z Eq v ). 



(e a + C) 
16tt 



(2\E 0y \ 2 ip + E* lz E 0y 



(7) 



Here Eo y is the y-component of the electric field am- 
plitude for the undistorted nematic (homeotropic orien- 
tation) and E\ z is the z-component of the field that is 
caused by nematic reorientation (calculated to the first 
order in ip). It is easily seen that in the undistorted LC 
the light maintains its polarization inside the layer, so 
that we have only one nonzero component of the elec- 
tric field E 0y (z) = E e lk * z , where k z — k^ e + ikj m ~ 
koy/s± — Sq + i-f±k n /(2y/e± — Sq) (terms of the order of 

(■~f±/{e± — Sq)) in k z are neglected because 7^ <C e±) 
and Eq is the amplitude of the incident electric field. In 
the linear approximation 9 remains zero. Straightforward 
calculations yield the following equation for E\ z (z) from 
Eqs. §|): ' 



(e_L +£a+ i{la + 7X)) 



d 2 E l2 



+ fco(e_L +«7-l) 



x(ej 



dz 2 

£a - Sq + z(7a + J±))E lz + k 2 (e a + «7a) 
d 2 (<pE Qy ) 



dz 2 



= 



(8) 



Substituting E\ z (z) into rt|) in the form Ei z (z) 



E(z)e 



and taking into account that koL 1 (L is 



the width of the layer), a first-order ODE for E(z) can be 
derived. Keeping in mind that E\ z (§) — we eventually 
obtain from Eqs. (0||) the following integro-differential 
equation for ip: 



dp(z,t) 



d 2 w(z,t) ( /7TK\ 







dt V 71 "/ 

^cos^— yz — z)J + sin ^— yz —z t 

< Z '-^V(2', t)dz' + ip(z, t)\ e - 2kl ™ z , 

where ip, £, k and r are parameters defined as: 

£a7-L - 3e a s±-f± + 2-y a e 2 L - 2-f a e±e a 



(9) 



2(e a + s±_)s a e± 
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e= 



27a£i - 3£ a £x7-L - £^7_l 



L 

K = — 



2 (e a + e±) £ e_L 



7^ 2 



7T 2^eT(e ± + e a ) ' tt 2 K 3 

In the parameters denned above only the linear terms in 
7a,7J_ were kept (70,71. <C 1). The parameter r is the 
characteristic time of the director motion and p = I/I c , 
where / is the intensity of the incident light and I c is 
defined as: 



Ic = 



IT 2 C(£_L +Sg)K 3 

L 2 Say/ell] 



V = (£a + 0/e a (10) 



J c coincides with the threshold intensity of the LIFT for 
a pure nematic (r] = 1, 71 = 711 = 0) at perpendicular 
incidence [^|. Then Eq. (^|) reduces to one obtained in 

We use a two-mode expansion with respect to z for the 
angle ip with the boundary conditions p(z = 0) = p(z — 
L) = 0: <p(z,t) = Ai(t)sin(7rz/L) + A 2 (t) sin (2nz/ L), 
where A\ and A 2 are time-dependent amplitudes. This 
is motivated by the fact that the distorted state is asym- 
metric with respect to the center of the layer because of 
absorption and the perturbation of the light polarization 
inside the layer. Therefore we have to include at least 
one mode that is symmetric and one mode that is anti- 
symmetric with respect to the center of the layer. After 
projecting Eq.(^) onto the trial functions we have a sys- 
tem of two equations for the modes A\ and A 2 : 



dA 1 dA 2 
t—j— = CuAi + L 12 A 2l t—j— 
at at 



C 21 A! +£ 22 A 2 ,{11) 



where the elements of the matrix £y depend on material 
parameters and the control parameters p and k (which is 
proportional to sjj). We look for solutions proportional 
to exp(at), where a is the growth rate. The procedure 
of deriving dj is straightforward but the expressions for 
these elements are too long to be presented here. 

The stability diagram in the (k, p) plane can now be 
calculated for any given material parameters of the LC. 
As an example we consider the nematic 5CB doped with 
the dye AD1 at 0.1% concentration. We used the fol- 
lowing values of material parameters at the temperature 
T = 24°: a a = 42 cm" 1 , n Q = 1.53, a e = 190 cm" 1 , 
n e = 1.71, (absorption coefficients and refractive in- 
dices of the ordinary and extraordinary light, respec- 
tively), A = 633 nm (wavelength of laser), £ = 58 ||, 
7 = 0.845 dyn ■ s/cm 2 , K x = 0.64 • 10" 6 dyn, K 2 = 
0.42 • 10~ 6 dyn, K 3 = 0.86 • 10" 6 dyn § ; the calcu- 
lations are made for a layer of 50 um thickness. For 
these parameters I c — 33.21 W/cm 2 , r = 2.49 s. It 
is easy to show the following relations: 7^ ~ a n /ko, 



7a — (a e n e -a o n o )/k , e± ~ n a , e a 



-rig (neglecting 



terms of the order (jj_/(sj_ — Sq)) ). 

The stability diagram is depicted in Fig. 2a. The solid 
line corresponds to a stationary bifurcation (Re(a) = 
Im(a) = 0) and the dashed one corresponds to a Hopf bi- 
furcation of the homeotropic state (Re(a) = 0, Ivnio ) 7^ 



0). These lines divide the (n,p) plane into a stable 
and an unstable region of the homeotropic alignment. 
They join in a so called Takens-Bogdanov point where 
det{C) = Tr{C) = 0. 

There are two differences compared to the case of a 
pure LC. First, the enhancement of the orientational op- 
tical nonlinearity described by the parameter £ leads to 
a " renormalization" of the threshold intensity (see Eq. 
(|l0|)). (However, since Fig. 2a is plotted with the renor- 
malized threshold intensity, this does not change the ap- 
pearance of the stability diagram). Second, absorption 
gives rise to the attenuation of the field inside the ne- 
matic. This results in a shift of the line of primary in- 
stability to the region of higher intensities as is shown 
in Fig. 2b. From this figure one can see the quantita- 
tive difference between the case when the absorption is 
neglected (dot-dashed lines) and when the absorption is 
taken into account (solid and dashed lines). Note that 
the critical intensity p t h for perpendicular incidence thus 
becomes larger than 1. 

It must be noted that we supposed that the nematic 
is maintained at constant temperature. Actually, due to 
the presence of the absorbing dye, the nematic will be 
heated by the light 0. We have estimated the max- 
imum temperature difference occurring inside the ne- 
matic 5CB doped with the dye AD1 from the steady- 
state heat conductivity equation (one-dimensional since 
we considered a plane wave). For the range of intensities 
/ = 30 — 100 W/cm 2 this difference was found to be no 
more than a few Kelvins. Thus we can usually neglect 
the temperature dependence of the material parameters 
and we took them to be constant across the layer. 



IV. STATIONARY DISTORTED STATE 

After the homeotropic state looses stability via a sta- 
tionary bifurcation (at not too large angle of incidence), 
the director settles in a stationary distorted state. One 
obtains the stationary z-dependent reorientation of the 
director by variation of the free energy density (jl]) with 
n defined in terms of the angles 9, tp. This gives us two 
ordinary second-oder differential equations for Q(z) and 
tp(z) (now these angles can be of arbitrary magnitude). 
We will not display them here because they are long and 
are obtained straightforwardly. 

These equations contain the field components which 
obey Maxwell's equations (^). It is convenient to write 
(^) using the Oldano formalism pT| . Following |l2| we 
write the matrix D (see (||)) as D = Dq + D z (z), where 





Vo 





1 

0/ 



(12) 



and e±_ = Sj_ + e a — e a + i"f a - The matrix D z 
contains the z-dependent angles 9 and tp. It is convenient 
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FIG. 2: a) Stability diagram of the homeotropic and station- 
ary distorted states in the (ft, p) plane. H is the region of 
the homeotropic state. SD is the region of the stationary dis- 
torted state bounded toward large p by the secondary Hopf 
bifurcation (dash-dotted line). TB is the Takens-Bogdanov 
point, b) Solid and dashed lines correspond to those in Fig. 
2a. The dot-dashed lines are obtained when the absorption 
effect is neglected, c) The secondary instability for small an- 
gles of incidence. Points A and B show the instabilities of the 
stationary distorted state for perpendicular incidence of the 
light with n x (z) — and n x (z) ^ correspondingly. 



to introduce a representation in terms of eigenfunctions 
of Do. The eigenvalue problem Doai = aion is solved by 
the eigenvalues 



a 2 = -ai = \ ej_ - sf. 



0,4 



-a 3 



(e_L + e a - sl)e±_ 
ej_ + e a 



(13) 



and eigenvectors 



ai,2 




«3,4 



/W e ± 
1 



V 



(14) 



We introduce the metric tensor 



M 



'0100 
10 
1 

,0 1 



(15) 



0, 



to define a scalar product between these vectors. With 
such a metric the eigenvectors are orthogonal to one an- 
other: cq Mctj = SijNi, where JV, is the "norm" of vector 
cii . The matrix Do is expressed by means of the vectors on 
as Do = Yli( a i/Ni)aiafM. The four vectors on give the 
polarization of four " proper" waves that propagate inside 
the layer without changing their state of polarization in 
the case of homeotropic alignment. The magnitudes of 
cii give the indices of refraction of these waves. Two of 
these vectors a% (a?) correspond to backward (forward) 
propagating ordinary waves and the other two 0:3 (04) 
correspond to backward (forward) propagating extraor- 
dinary waves. The contribution of the backward waves is 
negligibly small because the dielectric properties of the 
nematic change little on the spatial scale of the wave- 
length @. Thus we can expand ^(z) as follows 

*(») = b 2 {z)e tkoa2Z a 2 + b A {z)e ikoaiZ a4, (16) 

and write (pf) in terms of the amplitudes b 2 (z) and 64(2): 



= % [P22{z)b 2 + b 4 e-^ a *- a ^P 24 (z)] 



IMA J.V2 L - ' 

^ = ^[p 4i ( z ) bi + b 2 e 



—ik (a i —a 2 )z 



P*i(z)) , 



(17) 



where Pkj(z) = o. r [MD z (z)aj are the matrix elements of 
D z between the eigenvectors: 
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P22 
P 44 



e a e_L cos(6>) 2 sin(y) 2 

a l( e -L + e a COS(0) 2 COs(<^) 2 ) ' 2 

e a [(a4sin(6*) 2 — sin(20) cos(<^)s )( 



Pa 



e±(ej_ + e a )(e_L + e Q cos(6>) 2 cos(^) 2 ) 



[04 sin(0) — 003(6*) cos(iy9)so] sin(y) cos(6')e a 
fl2(e_L + e a cos(0) 2 cos(<£>) 2 ) 
e a )a 4 + e^s 2 (cos(6>) 2 cos(ip) 2 - 1)] 



(18) 



The advantage of the system ( p7| ) is that we now have 
only two equations for the "slow" amplitudes 62(2) and 
6 4 (z). So, we have a system of coupled ordinary differen- 
tial equations for 9{z), <p(z), 62(2) & ncl 64(2) with bound- 
ary conditions 6\ z =o,L = v\z=o,L = 0, and initial condi- 
tions b 2 \ z =o = A , &4| z= o = 0. Here A can be related to 
the normalized intensity p defined in the previous section: 
A Q = ^8tt 3 (s a + e±) (e± - s% + ij±) p/{e a s±ri). 

The system of "nematic+field" equations (with bound- 
ary conditions) is invariant under the transformation 
[6, tp, E x , E y ] — ► [8, —ip, E x , —E y ] owing to the reflection 
symmetry with respect to the y direction. Since the pri- 
mary instability breaks this symmetry, two different dis- 
torted states exist, which are mutual images under this 
transformation. For perpendicular incidence of the light 
there is an additional reflection symmetry with respect 
to the x direction and, as a consequence, the system 
of equations is also invariant under the transformation 
[0,<p,E x ,E y ]->[-0,<p,-E a ,E y ]. 

The system of equations can only be solved numeri- 
cally. For this purpose we introduced the new variables 
dd / dz, dip/dz to transform our set of equations to a sys- 
tem of six first-oder equations which was solved by the 
shooting method. To guarantee that we obtain the so- 
lution which originates from the homeotropic state we 
started with intensities only slightly above the thresh- 
old. Then, we increased p slightly and used the values 
of d8/dz\ z= ( h dip/dz\ z= Q obtained in the previous step as 
an initial guess. This procedure allowed us to derive the 
profiles 8(z), (p(z), 62(2) and 64(2) for any k and p above 
threshold. 

The director and field distributions for p = 2.0 and 
[3q = 11° (ft = 0.375) are shown in Figs. 3,4. In the next 
section we will show that for k — 0.375 the stationary 



distorted state becomes unstable at p c = 2.01, thus these 
figures represent the state slightly below the secondary 
instability. 



V. STABILITY ANALYSIS OF THE 
STATIONARY DISTORTED STATE 

Next we have performed a linear stability analysis of 
the distorted stationary state with respect to spatially 
periodic perturbations in the plane of the nematic layer. 
We write 

n = n (z) + Sn(x, y, z, t) = 
no(*) + 5n{z)e at+l ^ x+py \ * = * + * 1 = 
2 (bk{z) + Sb k (z)e at+i ^ x+py ^e ikoa " z a k , (19) 

k=2,4 

where Sri and Sbk are small spatially periodic perturba- 
tion with wavenumbers q and p; a is the growth rate. 

From the equation n 2 = 1 follows that n i5n = 0. 
Thus there are only two independent components of <5n. 
We obtained two linear equations for 8n x (z) and Sn y (z) 
which contain Sn x (z), 5n y (z) itself, their z derivatives up 
to second oder and 562,4(2) with complicated coefficients 
depending on the stationary distorted state 110(2), 62,4(2). 
Also, we decomposed the matrix D (see (||)) as D = Do + 
D z (z) + Di(5n), where the matrices Do, D z (z) correspond 
to the stationary state and were defined in the previous 
section, and the matrix Di(<5n) depends linearly on Sn. 
After linearization of (^) the equations for £62,4 can be 
obtained: 



d{Sb ^ z)) = %!■( Sb 2 P 22 + 5b 4 e iko ^- a ^ z P 2 4 + 6 2 P 9 ( , 1) + 6 4 e i M«4-a 2 ) 2! pa) 



N 2 



2-^22 



d(< y z)) = & ( Sb 2 e tk ^ a2 ~ a ^ z P 24 + 564P44 + b 2 e lka{a2 ~ a ^ z P^ > + b±P L 



24 
(1)' 



(20) 



where PL = a^MDiOj are the matrix elements of Di 
with respect to the eigenvectors ( |l4| ) and the Pjy were 
defined in ©. 

We have linearized the equation (^) substituting ip in 
the form (|19|). In principle we should have started from 



Maxwell's equations because the field perturbations con- 
tain x,y dependence. However this is a very good ap- 
proximation because the corrections are of the oder q/ko, 
p/k « 1. 

To solve the eigenvalue problem for a we expand 
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FIG. 4: Distortion of the field components inside the nematic 
layer for the stationary distorted state (p = 2.0, pa — 11°). 
E z is small compared to E x , E y and is not depicted; Eo is the 
amplitude of the incident electric field. 



^n-x,y(z), 5b2^{z) with respect to z in systems of func- 
tions which satisfy the boundary conditions (Galerkin 
method). For 5n the boundary conditions are 

8n x<y \ z —o,L = 0, thus we write Sn=^2 k Aksmfakz/L). 
Clearly the boundary conditions for the perturbations of 
the field amplitudes are <5&2.4|z=o = 0. One can see that 
at z = 0,L the r.h.s. of the system ( pfj| ) vanishes so 
one also has d(5b2,4)/dz\ z= o : L = 0. Therefore we used 
the expansion Sh = ^2 n B„ sin 2 (7rnz/(2L)). This set of 
functions is complete but not orthogonal. We have to 
truncate these expansions to a finite number of modes. 

We have solved the eigenvalue problem numerically to 
find the neutral surface p^{q,p) (for given angle (3q) which 
is defined by the condition Re(a(q,p)) = 0. The num- 
ber of Galerkin modes was chosen such that the accu- 



racy of the calculated eigenvalues was better than 0.1% 
(we took six modes for 6n and forty modes for Sh). 
The minimum of this surface gives the critical intensity 
p c =min 9iP po(q,p) and the critical wavevector (q c ,p c )- 
Since O c = Im(a) turned out to be nonzero at the min- 
imum, the instability corresponds to a Hopf bifurcation. 
The branch of the secondary Hopf instability is depicted 
as the dash-dotted line in Fig. 2a and for small angles 
of incidence in Fig. 2c. It is interesting to note the fol- 
lowing tendencies: as the incident angle (3q increases the 
critical intensity also increases, but the director and field 
deformations at the secondary instability decrease. 

The dimensionless Hopf frequency Q c t (r is defined in 
section III) versus k is shown in Fig. 5. Figure 6 shows a 
typical contour plot of the neutral surface po(q,p). The 
point (q c L,p c L) in this figure is the minimum of the sur- 
face and as is seen the bifurcation is inhomogeneous with 
some critical vector (q c ,Pc) ^ 0. This means that trav- 
elling waves are expected to appear. p c is only slightly 
below the homogeneous threshold po(q = 0, p — 0), which 
was calculated before for the pure LC [13) . 

As was pointed out in section IV, for nonzero (3q there 
are two symmetry-degenerate stationary distorted states. 
Clearly the two neutral surfaces are related by changing 
p to — p and the critical wave vectors will be (q c ,Pc) and 
(q c , —Pc)- Thus two different travelling waves with criti- 
cal vectors (g c ,±p c ) can be realized depending on which 
stationary state will be selected after the homeotropic 
state loses stability. 

An interesting situation arises in the limit of normal in- 
cidence. One might expect that for /3q — > the wavenum- 
ber q c — ► 0, since in this limit the external symme- 
try breaking in the x direction vanishes. However, this 
turned out not to be the case. The reason is that then 
another stationary instability that spontaneously breaks 
x-reflection symmetry intervenes the primary and the 
Hopf bifurcation. For the parameters of our computa- 
tion one has pth — 1-11, pel = 1.13 (point A in Fig. 
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FIG. 6: Contour lines for the surface p(p, q) correspond to 
fa = 11° (k = 0.375). The critical intensity is p c = 2.01 
with the critical wavevector (q c L,p c L) = (0.11, —0.06); po(q = 
0,p = 0) -p c = 1.5 -10" 3 . 



2c) and p c i = 1.17 (point B in Fig. 2c). One now has 
four symmetry-degenerate states and consequently four 
traveling waves with critical wave vectors (±g c , ±p c ). In 
some further investigations we have changed the ratios 
between the elastic constants keeping other material pa- 
rameters constant and saw the following tendency: the 
larger the anisotropy of the constants, the deeper the 
minimum of the surface becomes and the larger the mag- 
nitudes of the critical wavenumbers (see Fig. 7a, b). The 
absolute error of the dimensionless critical wavenumbers 
q c L,p c L depicted in this figure is less than 10~ 2 . In the 
one-constant approximation the bifurcation is homoge- 
neous {q c ,Pc — 0) for any n. This latter can be easily 
proved analytically. Perturbation theory can be used to 
investigate q, p dependence of the critical eigenvalue of an 
arbitrary stationary state. The calculations shows, that 
the perturbation of the eigenvalue is ~ p 2 + q 2 . 

From Fig. 7a, b one can see that q c L,p c L ~ 0.1. This 
means that the period of the structure 2ir /q c , 2it /p c ~ 
60L = 0.3 cm. Thus in an experiment the spot size of the 
light must be rather large in oder to observe the travelling 
waves. 

Finally we remark on the behaviour of the system in 
the nonlinear regime above the Hopf bifurcation. This 
system without transverse degrees of freedom has been 
studied extensively, and various regimes of complex be- 
haviour have been discovered. The bifurcation studied in 
this work marks the transition to simple periodic oscil- 
lations in the system without transverse degrees of free- 
dom, which is the first step towards complex behaviour. 
In models jl2| |l3| and simulations 0, a gluing bifur- 
cation was found above the secondary Hopf instability, 
which is a homoclinic bifurcation that restores the sym- 
metry broken by the Freedericksz transition. This glu- 
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FIG. 7: Critical wavenumbers q c ,Pc versus K2/K3 for differ- 
ent ratios K1/K3 (f3 = 11°). 



ing bifurcation was recently observed experimentaly [ jl5| . 
After this first gluing, complex nonlinear behaviour and 
eventually chaos was observed in both theory, simulation 
and experiment []l6|| . An analogous gluing bifurcation 
should exist also in the case of the spatially extended 
system. 

The behaviour of the system in the vicinity of this glu- 
ing bifurcation, can, however be radically different from 
what was observed in the experiment |J15| . In the spa- 
tially constrained system (i.e. the director oscillation 
induced by a narrow beam as observed in the experi- 
ments) one observes stochastic behaviour in the vicinity 
of the first gluing only as a consequence of experimental 
noise. It has been shown, however, Jl7| that any spa- 
tially extended system, which possesses a homogeneous 
limit cycle (which is stable with respect to homogeneous 
perturbations) becomes unstable as it approaches a ho- 
moclinic bifurcation. This instability is either a phase 
instability, or a finite-wavelength period-doubling insta- 
bility. On these grounds one can expect to observe very 
complicated behaviour (probably spatio-temporal chaos) 
in our system already at the threshold intensity of the 
first gluing. As opposed to the previous case, this would 



8 



be true deterministic chaos, not merely stochasticity due 
to noise. 



VI. CONCLUSION 

We have found the threshold of the LIFT for the 
homeotropic state and the threshold of the secondary in- 
stability of the stationary distorted state in a nematic 
LC, including the dye-doped case, for different incident 
angles of the light. In particular we have demonstrated 
that the stationary distorted state loses stability in an 
inhomogeneous Hopf bifurcation with some nonzero crit- 
ical wavenumber that leads to the formation of travelling 
waves in the plane of the layer. 

Our result demonstrates a general feature of Hopf bi- 



furcations in spatially extended systems with broken re- 
flection symmetry, as is the case in the LIFT-distorted 
state. Except for special cases, like those where the re- 
flection symmetry can be restored by going into a moving 
frame, the neutral surface exhibits the signature of the 
broken symmetry. Consequently, at q c = p c = the neu- 
tral surface does not have a stationary point and thus 
cannot have a minimum. This general feature was ap- 
parently first noted in the context of reaction-diffusion 
systems JL8| . 
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